Measuring gravito-magnetic effects by multi ring-laser gyroscope 
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We propose an under-ground experiment to detect the general relativistic effects due to the 
curvature of space-time around the Earth (de Sitter effect) and to rotation of the planet (dragging 
of the inertial frames or Lense-Thirring effect). It is based on the comparison between the lERS value 
of the Earth rotation vector and corresponding measurements obtained by a tri-axial laser detector 
of rotation. The proposed detector consists of six large ring-lasers arranged along three orthogonal 
axes. In about two years of data taking, the 1% sensitivity required for the measurement of the Lense- 
Thirring drag can be reached with square rings of 6 m side, assuming a shot noise limited sensitivity 
(20prad/s/\/ Hz). The multi-gyros system, composed of rings whose planes are perpendicular to 
one or the other of three orthogonal axes, can be built in several ways. Here, we consider cubic 
and octahedron structures. The symmetries of the proposed configurations provide mathematical 
relations that can be used to study the stability of the scale factors, the relative orientations or the 
ring-laser planes, very important to get rid of systematics in long-term measurements, which are 
required in order to determine the relativistic effects. 

PACS numbers: 42.15.Dp, 42.30.Sy, 42.55.Lt, 91.10.Nj 
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I. INTRODUCTION 

The general theory of relativity is the most satisfactory 
description of gravitational phenomena. The theoretical 
breakthrough came with Einstein's geometrical represen- 
tation of gravity: as different test masses fall in the same 
way in a gravitational field, gravity must be a property 
of space and time rather than of the masses themselves. 

Until now, almost all successful tests of general rel- 
ativity (Shapiro time delay [5^, light deflection by the 
sun [6], perihelion shift of the orbit of Mercury [7]) have 
been probing the gravitational field of the Sim, without 
considering its proper rotation. However, general relativ- 
ity predicts that the stationary field of a rotating body 
is different from the static field produced by the same 
non- rotating mass. The difference is known as gravito- 
magnetism and consists of a drag of space-time due to 
the mass currents. The rotational frame-dragging effect 
is also known as the Lense-Thirring (LT) effect. 

A direct experimental evidence of the existence of the 
GM field has been obtained so far by Ciufolini [9] and 
by Francis Everitt and the GP-B group [3]. The Lense- 
Thirring effect, averaged over several orbits, has been 
recently verified by analysing the node orbital motion of 
two laser ranged freely falling satellites (LAGEOS-1 and 
LAGEOS-2) which orbit the Earth. In the measurement 
presented in Ref. the two LAGEOS satellites were 
used to confirm the LT effect with an accuracy of the 
order of 10%. However, the launch of a third properly 
designed satellite LARES will give the opportunity to 
measure the LT effect with an accuracy of the order of 
1% (citare). 

The possibility to detect Lense-Thirring with ring 
lasers has been discussed in the past [Ml [53] . Recently it 
has been already pointed out that a multi-gyros system 
is able to test locally the Lense-Thirring effect [10]: an 
array of six, 6 m side, square ring-lasers have enough sen- 
sitivity for this purpose. The rings must have different 
orientation in space. In the present paper we concentrate 
the attention on the symmetries of the rings arranged 
on the faces of a cube or along the edges of an octa- 
hedron, extracting the relevant relations important for 
the diagnostics of the system. At the end we summarize 
and sketch the proposed experiment. For completeness 
we must mention that an experiment of the type we are 
planning and preparing could also be made in principle 
using matter waves instead of light. This possibility has 
been proved experimentally for various types of parti- 
cles such as electrons [T2], neutrons [13], Cooper pairs 
[Hj . Calcium atoms [15], superfiuid He3 [TH and super- 
fluid He4 [T7] . Cold atoms interferometry, in particular, 
yields very high sensitivity and it is suitable for space 
experiments because of the apparatus small size. How- 
ever, atoms interferometry experiments in space do not 
provide an independent measurement of the Earth an- 
gular velocity, are affected by the mass distribution of 
the Earth, and test the average of the relativistic effect 
rather than the local one. Eventually, the comparison 



between in-space and on-ground measurements could be 
very valuable. 

II. DETECTION OF GRAVITO-MAGNETIC 
EFFECTS 

Gravito-magnetism (GM) is a general relativistic phe- 
nomenon related to the presence of mass currents in the 
reference frame of a given observer. In the case of celestial 
bodies, including the Earth, and excluding translational 
motion with respect to the center of the body, gravito- 
magnetic effects are due to the absolute rotation of the 
massive source with respect to distant stars. When the 
Einstein equations in vacuum are applied to this kind of 
symmetry and are linearised (weak field approximation) 
GM is accounted for by the analogue of a magnetic field 
of a rotating spherical charge. In practice at the lowest 
approximation level, a dipolar GM field is obtained, with 
the dimensions of an angular velocity. Its explicit form 
in a non-rotating reference frame centred on the source 
(in our case the Earth center), is (see e.g. [46] ) 

2G 

^ ^ t"^® ~ ^('^ffi ■ '^r)Ur] (1) 

where R = Rur is the position of the laboratory with 
respect to the center of the Earth and J® is the angu- 
lar momentum of the Earth, whose modulus is of course 
given by the product of the moment of inertia of the 
planet multiplied by its angular velocity. 

The effect produced by a field like ([!]) on a massive test 
body moving with velocity v looks like the one produced 
by a magnetic field on a moving charge: in fact, the 
geodesic equation in weak field approximation reads 

'l-=G + vAB (2) 

where G = ~GM / R^Ur is the Newtonian gravitational 
field, so that the effect can be described in terms of a 
gravito-electromagnetic Lorentz force, where the New- 
tonian gravitational field plays the role of the gravito- 
electric field (GE). 

Furthermore, the rotation of the source of the gravita- 
tional field affects a gyroscope orbiting around it, in such 
a way that it undergoes the so-called Lense-Thirring pre- 
cession, or dragging of the inertial frames of which the 
gyroscope defines an axis [151 IMl [55] . This phenomenon 
shows up also when one considers a freely falling body 
with local zero angular momentum (ZAMO: Zero Angu- 
lar Momentum Observer): it will be seen as rotating by 
a distant observer at rest with the fixed stars [TS] . 

A. Mechanical gyroscopes 

Gravito-magnetic effects can in principle be measured 
applying different methodologies. The one that has most 
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often been considered is focused on the behaviour of a 
gyroscope, that can be either in free fall (on board an 
orbiting satellite) or attached to the rotating Earth. The 
axis of the gyroscope is affected in various ways by the 
presence of a gravitational field. As for GM, a little me- 
chanical gyroscope is the analogous of a small dipolar 
magnet (a current loop), so that it behaves as magnetic 
dipoles do when immersed in an external magnetic field. 

When studying the motion around the Earth of a gy- 
roscope whose spin vector is S , one is led to the formula 

[zms]: 



dS 



= n' AS 



(3) 



In Appendix [A] we work out the explicit expression of ft' 
in general relativity and, more in general, in metric the- 
ories of gravity, using the Parametrized Post-Newtonian 
(PPN) formalism [48] : we show that it is related to the 
gravito-magnetic components goi of the metric tensor 
and its expression is given by (see Eqs. (A6)-(A10)) 
ft' = flc + rJs + riw + ^T, so that we can distin- 
guish four contributions, namely the geodetic term ilc, 
the Lense-Thirring term fls, the preferred frame term 
Jln/, the Thomas term Jly. All terms in H' are called 
relativistic precessions, but properly speaking only the 
second is due to the intrinsic gravito-magnetic field of 
the Earth, namely it is CIb — ~\B, and manifests the 
Lense-Thirring drag. 

Ciufolini j2| deduced the relativistic precession of the 
whole orbital momentum of two LAGEOS satellites 
whose plane of the orbit is dragged along by the rotating 
Earth. Again on Eq. (§ was based the GP-B exper- 
iment, whose core were four freely falling spherical gy- 
roscopes carried by a satellite in polar orbit around the 
Earth I? . While time goes on and the available data grow 
it is expected that the Lense-Thirring drag will emerge 
from the behaviour of the unique (so far) double pulsar 
system [4]. 



B. Using light as a probe 

A different experimental approach consists in using 
light as a probe. In this case the main remark is that 
the propagation of light in the gravitational field of a 
rotating body is not symmetric. The coordinated time 
duration for a given space trajectory in the same sense as 
the rotation of the central source is different from the one 
obtained when moving in the opposite direction. This 
asymmetry would for instance be visible in the Shapiro 
time delay of electromagnetic signals passing by the Sun 
(or Jupiter) on opposite sides of the rotation axis of the 
star (or the planet) pUj^T]. 

This property of the propagation of light is the one 
which we wish to exploit in our Earth-bound experiment 
using a set of ring lasers. In a terrestrial laboratory, light 
circulating inside a laser cavity in opposite directions is 
forced, using mirrors, to move along a closed path in 



space. What is closed from the view point of the labo- 
ratory is not so for a fixed-stars-bound observer, but the 
essential is that the two directions are not equivalent and 
that the two times required for light to come back to the 
active region are (slightly) different. As it happened al- 
ready in the case of the mechanical gyroscopes, here too 
the difference in the two times of flight is made up of var- 
ious contributions depending on the rotation of the axes 
of the local reference frame with respect to distant stars, 
on the fact that the local gravitational (Newtonian) po- 
tential is not null, and of course on the GM drag (which 
is our main interest). What matters, however, is that the 
final proper time difference (a scalar quantity) is invari- 
ant: it does not depend on the choice of the reference 
frame or of the coordinates. 

Performing the calculation in linear approximation for 
an instrument with its normal contained in the local 
meridian plane (see the Appendix |A] details) we find 



c5t 



— "a 
c 



2— — smfc'sma 



(2 cos 9 cos a -t- sin 9 sin a) 



(4) 



where A is the area encircled by the light beams, a is the 
angle between the local radial direction and the normal 
to the plane of the instrument, measured in the meridian 
plane, and 9 is the colatitudc of the laboratory; is 
the rotation rate of the Earth as measured in the local 
reference frame (which includes the local gravitational 
time delay). 

Eq. Q can also be written in terms of the flux of an 
effective angular velocity VI through the cross section of 
the apparatus: 



5t=-A- O, 



(5) 



where A = Aun is the area enclosed by the beams and 
oriented according to its normal vector m„. In particular, 
it is $7 = -|- f2', and the term proportional to Jl® is 
the purely kinematic Sagnac term, due to the rotation of 
the Earth, while fJ' — Q,g + + + encodes the 
relativistic effects (see Appendix [A]) 

For a ring laser in an Earth-bound laboratory, the 
geodetic and Lense-Thirring terms are both of order ~ 
10~^ with respect to the Sagnac term, while the Thomas 
term is 3 orders of magnitude smaller. As for the pre- 
ferred frame term, the best estimates |56l ETj show that 
this effect is about 2 orders of magnitude smaller than 
the geodetic and Lense-Thirring terms. Consequently, to 
leading order, the relativistic contribution to the rotation 
measured by the ring laser turns out to be $7' ~ rJc+ris, 
which we aimed at measuring in our experiment. In other 
words, the goal of our experiment will be the estimate of 
Jl' (see Fig. [ij which embodies the gravito-magnetic ef- 
fects in a terrestrial laboratory. 

In particular, the proposed experiment can also pro- 
vide high precision tests of metric theories of gravity 
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FIG. 1: The amplitude of the relativistic effects on the sur- 
face of the Earth, according to the theory of general relativity, 
in units oiprad/s, as a function of the colatitude 6. The con- 
tinuous, dashed and dotted lines correspond to f2' = flc + fis 
projected along the directions: i) parallel to fJ® (i.e. ii) 
u-r (local radial or zenithal direction); and iii) ue (local North- 
South direction), respectively. To evaluate the contribution 
to r2' from fJc, we have projected ilc along fJ® (continuous 
line plus triangles) and ug (dotted lines plus squares). We 
note that the gravito-electric term has only the ue compo- 
nent and therefore along the radial direction we have a pure 
gravito-magnetic term. 



which are described in the framework of (PPN) formal- 
ism. In fact, from Eqs. ( |A22|A23[ ), we see that, on set- 
ting for the rotating Earth J — 70^20, we obtain 



VIg = -(1 + 7)^^ sini9rj0Mfl, 



nn = - 



,GM 

4 



1+7+^ G/ff 



(6) 

3 (r20 • Ur) Ur] (7) 



where ai and 7 are PPN parameters (e.g. ai — Q 
and 7 = 1 in general relativity) which account for the ef- 
fect of preferred reference frame and the amount of space 
curvature produced by a unit rest mass, respectively. 

As shown in Sect. |VIIB from a high precision mea- 
surement of the vector H' in the meridian plane, we 
should be able to place new constraints on the PPN pa- 
rameters ai and 7. 



III. THEORY OF THE MEASUREMENT: 
COMBINING TOGETHER THE RESPONSE OF 
SEVERAL RINGS 

A. The response of a ring laser 

A ring laser converts time differences into frequency 
differences. In fact, since the emission is continuous, the 
right handed beam adjusts itself to give a standing wave 
whose wavelength is an integer sub-multiple of the space 



length of the loop P: ct+ = P = iVA+. The same hap- 
pens with the left handed beam, but being the total time 
different, also the wavelength of the corresponding stand- 
ing wave will be different: ct„ = 7VA_. The two modes 
of the ring can have different N, a situation usually called 
'split mode', but the higher accuracy of the measurement 
has been obtained so far with the two modes with equal 
N. Considering the time of flight difference in terms of 
the wavelengths of the two standing waves we see that: 



c5t = iV(A+ - A_) = Nc 



f--h 



P 



PA 



51 



(8) 



The ring laser equation [22] relates the frequency split- 
ting 5f of the two optical beams inside the ring interfer- 
ometer with the experienced rotation rate of its mirrors 



(9) 



where P is the perimeter and A is the laser wavelength. 
The response i? of a ring laser to the rotation rate Jl, in 
units of rad/sec, is simply a rescaling of the frequency 
splitting by the scale factor S* = 44, i.e. 



R = Sf/s = u„ • n. 



(10) 



The scale factor S plays a crucial role in the accuracy 
of the measurement of ft and to estimate the relativistic 
effects the ratio |p must be known and kept at 10^^" 
accuracy level for months. The requirements to keep the 
apparatus in the optimal working conditions will be dis- 
cussed in section Hvl 

Since the effective angular velocity as well as the 
gravito-magnetic one is of the order of lO~^n0, angles 
between vectors must be measured at the corresponding 
accuracy level. Unfortunately, the absolute measurement 
of Un in the fixed stars reference system with the accu- 
racy of nano-radians can hardly be achieved. However, 
we can relax this requirement by using M > 3 ring lasers 
oriented along directions u" (a = 1 . . . M), where not all 
lie in the same plane. In fact, Jl can be completely 
measured by means of its projections on at least 3 inde- 
pendent directions (e.g. defining a tri-dimensional Carte- 
sian system) and the redundancies of the measurement 
can be used as a monitor and control of the stability of 
directions w". We further assume that ring lasers have 
identical sensitivity and noise parameters. From an ex- 
perimental point of view this can be easily satisfied by 
building the devices with scale factors that differ less than 
%. 

In order to simplify the sensitivity calculations of the 
system one can consider multi-axial configurations en- 
dowed with symmetries. As all the ring laser normals u" 
are equivalent in space, symmetric configurations should 
be more efficient in the rejection of spurious effects and 
in the control and monitoring of the relative orientation 
of the ring lasers. The natural choice is to take advan- 
tage of space symmetries of regular polyhedra, setting 
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one ring for each plane parallel to their faces. If we do 
not consider the degeneration between opposite faces, we 
have M=3 in the case of the cube, 4 for tetrahedron and 
octahedron, 6 for dodecahedron, and 10 for icosahedron. 
There is a peculiar geometry with M = 3, obtained by 
arranging the rings along the edges of an octahedron, 
where the different rings can be nested together, sharing 
2 by 2 the same mirrors. We will refer to it in the fol- 
lowing by speaking of "octahedral configuration". The 
M = 3 is the minimum number of rings necessary to re- 
construct the rotational vector, but a redundancy is very 
appropriate to enhance statistic and to have control tests 
on the geometric accuracy. 

In general, by simple arguments, one can demonstrate 
that, for regular polyhedra configuration 



M 



^ tt„ = {M > 3) 



and that 



M 



M 



(M> 3). 



(11) 



(12) 



As a consequence, one can study linear and quadratic 
combinations of ring lasers responses Ra which are invari- 
ant under permutations of the ring laser labels a, i.e. L = 
Ra and Q = J2a ^a- -f'or non-symmetric configura- 
tions we can generalize their definition as L = LaRa 
and Q — QapRaRp , where La and Qap are suitable 
constants which depend on Ua- The interest in such lin- 
ear or quadratic forms relies on their behaviours in the 
presence of noise fluctuations or variations of the geom- 
etry of the configuration. 

They also allow us to carry out analytical estimates of 
the overall sensitivity of a tri-axial system of ring laser 
to relativistic effective rotation rates. 



B. Requirements for the geometry of the 
configuration 

The response of each ring laser can be conveniently 
written as 



Ra^^ 



SUa) 



(13) 



where Sua = SSa Ua + Si9a A Ua account for system- 
atic errors in the scale factors and orientations in space 
and Ea represents the additive noise that affects the ro- 
tation measurement Ra , that we assume averaged on the 
observation time T ~ 1 day. We assume as well that 
Ea are Gaussian distributed random variables with zero 
mean and variance ctq. Modulus |<5itQ.| ~ (55*0 and di- 
rection di9a A Ua represent the deviations from regular 
polygon geometry in the plane and from polyhedra ge- 
ometry in the space, due to scale factor fluctuations 6Sa 
and infinitesimal rotations S'da, respectively. In what 



follows the crucial assumption is that systematic errors 
(scale factors and relative alignment of Mq) are negligi- 
ble with respect to statistical errors, i.e. \fl ■ Sua\ < ctq 
or equivalently \Sua\ < aci/fl, while the dihedral angles 
arccos(itQ, -it^) can nearly approximate a regular polyhe- 
dron configuration. 

Redundancy of responses, if M > 3 rings are involved, 
can be used to control systematic errors projected along 
the direction oi CI. In fact, the rigidity of the configu- 
ration imposes some linear kinematic constraints among 
different estimates of the laboratory rotation. In general, 
any linear combination of 3 responses Ra gives an esti- 
mate of the local rotation fl and we can test the consis- 
tency among different estimates by means of the ordinary 
least square fit. A very simple linear constraint can be 
found for regular polyhedral configurations 



M 



L — ^ Ra 



(14) 



and we will illustrate its statistical property as an exam- 
ple of the power of the method. 

From the definition of L immediately follows that it is 
Gaussian distributed with zero mean and standard devi- 
ation (Tl ~ y/Man- In addition, possible misalignments 
di9a/\Ua or scale factor fiuctuations SSaUa are amplified 
by a factor of in the mean value of L 

M 

<L> = ^n-Sua (15) 

M 

= nY{5SaUa+&'0af\Uaf , (16) 



without affecting the corresponding variance a\. Thus 
< L > can be used as a "null constraint" which is mini- 
mum when the configuration geometry is a regular poly- 
hedron, and so the overall mean error parallel to ft can 
be monitored at ^ ^/Man/^ — 10^^" accuracy level. 



C. Estimate of the parallel component of the 
relativistic effective rotation vector 

An estimate of fl^ for symmetric configurations readily 
follows from Eq. ( 12 ) 



M 



(17) 



a=l 



r. M M 

- ' E-""-" + l^E^^^ (18) 



M 



Q — 1 Ct — 1 

Its mean value and standard deviation read (see App. [B 

1 . 



< g > 



18 4 12^, p 



(19) 
(20) 
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In addition, one can demonstrate that Q is non-central 
distributed with M degrees of freedom and non-centrahty 
parameter fi^. In order to estimate the relativistic effec- 
tive rotation, we must subtract fi® from the rotation rate 
estimated in the laboratory. To this end we calculate the 
difference A ee {Q ~ n%) that, in the limit of high SNR 
(|r2|/crn >> 1)7 tends to be Gaussian distributed with 
mean 



and standard deviation 



CIA ^ (2V3/\/m) 17e an 



(21) 



(22) 



where we have neglected terms of the order of ao/il. 
The SNR =< A > /cta of the parallel component of 
relativistic effective rotation is increased by a factor of 
^JmJZ with respect to the sensitivity of each ring laser. 

The advantage of this approach is that we compare 
scalar quantities (moduli of rotation vectors) measured 
with respect to the local and distant stars reference sys- 
tems. Its drawback is the very poor sensitivity to the 
perpendicular component ilj^ of the relativistic effective 
rotation. In fact, fi^ - fi^ = 211' • Jl® -l-|r2'P, and the 
ratio between the second term (which is associated to the 
perpendicular component as |r2'p = il'j^ + O^) and the 

first term is ~ GM/c^R ~ IQ'^^ . 

It is worth noticing that statistical fluctuations of L 
(control of geometry by redundancy) and Q (measure of 
relativistic effects) are uncorrelated, and that they tend 
to be independent in the limit of high SNR. 



D. Estimate of the components of the relativistic 
effective rotation vector 

By arranging the response of ring lasers Ra as 
M-tuples in a M-dimension vector space R = 
(i?i, i?2, ^3, • • • 7 Rm), we can easily define projection op- 
erators that allows the estimate of local meridian plane 
Ai and also the direction to of in the physical space. 
Moreover, the norm of projected random vectors are de- 
scribed by remarkably simple statistics. According to the 
definition of the matrix product we have R = TV $7 -I- e, 
where A/' is a M x 3 matrix whose elements are N^i — 
{ua)i and e = (ei, £2, £37 • • • 7 £m)- Thus, the random vec- 
tors R can be projected on the linear subspaces Vm and 
Qm of dimensions 2 and M — 2, which represent respec- 
tively a plane in the physical space and its complemen- 
tary space. The physical symmetry of the rotating Earth 
imposes that the relativistic effective rotation vectors and 
fl^ lie in the same plane, i.e. the meridian plane, and 
therefore the knowledge of the orientation of this plane is 
crucial if we want to measure not only the modulus but 
the whole vector. We recall that a plane is defined as the 
set of the points sv 4- tw, where s and t range over all 
real numbers, v and w are given orthogonal unit vectors 
in the plane. 



The parallelism of v A w with the normal to the 
meridian plane can be tested under the hypothesis that 
the rotation signal is fully located in the Vm subspace 
while the Qm subspace contains only noise. The test 
can be easily performed over the norms of the two 
projections Ep{v,w) = and Eq{v,w) = 

I |(5„,.u,il| where we have introduced the symbol | |i?|| = 

(Eaii^a)^^^ to indicate the L-2 norm in the M- 
dimensional Euclidean response space. The M x M pro- 
jection matrices P^^w and Qv,u} can be written explicitly 
as functions of the unit vectors v and w 



Qr 



NV(NVf 
I -NV{NVf 



(23) 
(24) 



where V is a 3 x 2 matrix with columns v and w, and 
/ is the M X M identity matrix. As shown in appendix 
[b] , the probability distribution of Ep is non-central 
with 2 degrees of freedom and non-centrality parameter 
f2^, while Eq should be distributed with M-2 degrees 
of freedom. 

The best estimate of v and w is then obtained by 



(u, w) 



arg max 

v.w 



\Pv.wR\ 



(25) 



where the taken over the unit sphere. The 

direction of the Earth rotation axis can be estimated as 



a particular case of Eq. (10 1. In fact, the projectors P^ 
and Qw can be obtained by substituting the matrix V 
for the 3x1 matrix W with columns w. The difference 
lies in the dimension of the corresponding subspaces, i.e. 
■pvv and Qyy have dimension 1 and M-1 respectively. It 
is worth noticing that the maximum of Eq. ( 25 1 can 



be computed by an analytical formula both for the loca- 
tion of the meridian plane and the direction of the Earth 
rotation axis. In fact, if we introduce in the local ref- 
erence frame the (local) spherical coordinate {R, 8 and 
$) (we use capital letters to avoid confusion with Sect. 
II) and parametrize the unit vectors v and w with these 
angles, for instance w = (cos $ sin O, sin $ sin 9, cos 8) 
and V — (cos $ cos 8, sin $ cos 8, — sin 8) , we have that 
the maximum of and ||Pu)-R|P is achieved for 



tan 8 = 
tan<l) = 



Whr) 
r^kr V'^ 

RTJR I 



(26) 



where F, H, K, J are M xM symmetric matrices which 
are functions of the Ua alone. 

In general, there are no analytical calculations for 
mean and variance of v, w and one must run Monte Carlo 
simulations to get their estimates. However, in the limit 
of high SNR Ep and Eq tend to be Gaussian distributed, 
as well as fluctuations of v and w around their mean val- 
ues. The same reasoning holds true also for the estima- 
tion of 8 and $. 
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The validity of the proposed experimental configura- 
tion has been checked by a numerical simulation over a 
period of 1 year of the six responses of the octahedral 
configuration oriented as in Fig. [23] of subsection |V C| 
In order to simplify the calculations we assume that the 
laboratory colatitude \s 9 — tt/A and that the normal 
to the plane of a ring forms a 7r/4 angle with respect 
to the Earth axis, and another normal is orthogonal to 
the former and forms again a 7r/4 angle with the west- 
east direction. This configuration is close to a possible 
experimental arrangement at the Gran Sasso National 
Laboratories (LNGS) within few degrees. The directions 
of the unit vector Ua in the local reference frame are 
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and the rotation signal for the 6 rings are equal within 
a factor -^2. We assume one mean sidereal day Ts = 




-0.4 -0.3 -0.2 -0.1 0.0 
rotation rate difference [prad/s] 



FIG. 2: Histograms of the difference A between Q and Jl^, 
normalized with the mean sidereal day, collected for 3 months 
(dark histogram) and one year (light histogram). 

86164.0989 s of integration and a noise standard devi- 
ation aQ — 7 X 10~^ prad/s. The variance ctq is ex- 
trapolated from present "G" sensitivity at 10'* s and 
scaling by a factor 5, due by the increase of the ring 
size and the power of the laser of a factor 1.5 and 
10, respectively. The relativistic rotation contributions 
n'^ = -2.8 X IQ-'^prad/s and n'g = -5.6 x IQ-'^prad/s 
have been added to the Earth rotation vector fl^ , as esti- 
mated by lERS [2S] . The component of relativistic effects 
parallel to is n\, = {n'g + n'^)/\/2 = 5.9xl0''^prad/s. 



Using Eq. ( 10 ) we calculated the responses of the 6 rings 
and then we injected the Gaussian noise. In Fig. [2] we 
show the histograms of Ts A/ {2Tr) accumulated for 90 
and 366 sidereal days. The corresponding mean val- 
ues of the parallel component of relativistic effects are 
—6.0 X 10~^ prad/s and —6.2 x 10~^ prad/s with stan- 
dard deviations 4.7 x 10^^ prad/s and 2.6 x 10"'^ prad/s, 



respectively. Thus a ^ 10% accuracy can be achieved 
in 3 months by simply comparing the square modulus 
of rotation vectors. In order to give a full estimate of 
the vector fl' , we have also explicitly calculated day by 
day the angles 8 and $ for describing the orientation of 
the meridian plane and the direction of the Earth rota- 
tion vector. The results are summarized in Fig. [3] and [4] 
where we report the time evolution of these angles and 
in Fig. [5] where we show the corresponding annual polar 
motion. 



2.5x10"^ 



^ 2.X10- 
to 



(27) ^ 
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Sidereal days 

FIG. 3: Change of the angle O due to polar motion as mea- 
sured by the ring laser responses in one year. 

By synchronizing the polar motion measured in the lo- 
cal reference system with the polar motion measured by 
lERS in the fixed star reference system, the two reference 
frames will coincide within the accuracy of the measure- 
ment of O and say 1 part of lO^''. As a final remark, 




50 100 150 200 250 300 350 

sidereal days 

FIG. 4: Change of the angle $ due to polar motion as mea- 
sured by the ring laser responses in one year. Note the large 
variation of $ which correspond to a nearly complete preces- 
sion cycle of the Earth axis in one year. 

we point out that the full measurement of the vector J7' 
allow us for the estimate of il^ ~ 2 x 10~^ prad/s with a 
standard deviation of the same order of magnitude of the 
estimate of J7J| . This represent an increase of the relativis- 
tic rotation signal of ~ 30 %. However, the estimate of 
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n' is crucial to separate the geodetic from Lense-Thirring 
contributions and/or to measure the PPN parameters ai 
and 7. 

0.5| 




0.1^ - 

-0.1 0.0 0.1 0.2 0.3 

yu rad 

FIG. 5: The estimated polar motion from the 6 ring laser 
responses. 



E. The Earth motion and feasibility of the 
experiment 

Since our goal is the estimate of the Lense-Thirring ef- 
fect at few % accuracy, the independent measurement of 
-|- ^jfEL ^ which represents the rotation of the labora- 
tory with respect to distant stars, must be determined to 
lO~^°r20. Due to tidal forces and to the exchange of an- 
gular momentum between the solid Earth and geophysi- 
cal fluids, the angular velocity of the Earth varies in time, 
both in direction and modulus. Changes in modulus cor- 
respond to a variation of the Length of the Day (LoD) 
of few milliseconds with respect to atomic clocks. The 
direction of the rotation axis of the Earth varies with re- 
spect to both the fixed stars and the Earth-fixed reference 
frames. Nowadays, the best Earth rotation monitoring is 
provided by the lERS 05C04 time series [28] which are 
routinely obtained using the geodetic space techniques 
VLBI (Very Long Baseline Interferometry) , SLR (Satel- 
lite Laser Ranging), GPS ( Global Positioning System) 
and DORIS (Doppler Orbitography and Radioposition- 
ing Integrated by Satellite). 

In Figs. |6] and [7] we report the Length of the Day 
(LoD) and the pole position with the corresponding er- 
rors of the last six years. It is worth to noticing that the 
achieved precision is 0.001 ms in the LoD and 0.1 ^arcsec 
in the pole position. 

Further improvements are expected in the next few 
years and the overall errors in LoD and pole position 
should decrease of a factor 10 that is crucial for a 1% 
measurement of the relativistic rotation terms. However, 
the lERS 05C04 time series is already sufficient to get 
jfl®! with 3% accuracy. 

For what concerns the differential rotation of the lab- 
oratory with respect to the rotation estimated by lERS, 
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FIG. 6: The change of the Length of the Day (LoD) over the 
last 6 years from the lERS 05C04 time series. Notice that 
estimated errors of LoD decreased in the last years to a level 
which correspond to 10~^* rad/ sec, i.e. 0.1 ppb £7^. 
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FIG. 7: The change of the direction of the Earth rotation 
axis (i.e. pole position) over the last 6 years from the lERS 
05C04 time series. Estimated errors are also plotted. 



it is expected to be sufficiently small to contribute to fi^ 
only trough ^(^\\ ■ However, rj^^j^ is still largely unknown 
due to possible micro-rotations of the crust of the Earth. 
This is one of the causes limiting the performances of G in 
Wettzell: the Earth crust motion caused by atmospheric 
changes. It is assumed that an underground facility is 
less sensitive to this kind of noise sources. It is as well 
important to keep the experiment close to VLBI stations. 
The underground Gran Sasso Laboratories is placed half 
way between two relatively close VLBI stations, Medic- 
ina and Matera [SD] which can provide estimates of 
the crustal motion of the Adriatic plate [31]. A signifi- 
cant contribution to O^^^ comes from the "diurnal polar 
motion" (periodic motion of the Earth crust due to tides) 
and consists in periodic changes of amplitude ^ lO^^fJ^. 
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This effect has been already measured by large ring laser 
gyroscopes [35] , and can be accurately modeled and then 
subtracted from ring lasers responses. 

We conclude that by means of available geodesies and 
geophysics techniques, provided that the experiment is 
located in an area with very low relative angular motion 
(^REt)j ^ suitable tri-axial detector of rotation can in 
principle detect il' with % precision. 



IV. THE 'REAL APPARATUS', THE PRESENT 
SENSITIVITY OF G IN WETTZELL 

Sensor properties 

A closer look at equation [9] reveals that there are three 
basic effects one has to carefully account for. These are: 

• scale factor stability {4A/XP) 

• orientation of the gyroscope with respect to the in- 
stantaneous axis of rotation of the Earth 

• instantaneous rate of rotation of the Earth - Length 
of Day (LoD) 

The scale factor for all practical purposes has to be held 
constant to much better than 1 part in 10^". Otherwise 
the frame-dragging parameter cannot be determined un- 
ambiguously. For G, the base of the gyroscope has been 
manufactured from Zerodur, a glass ceramic with a ther- 
mal expansion coefficient of a < 5 x 10~^/°C. Further- 
more the instrument is located in a thermally insulated 
and sealed environment with typical temperature varia- 
tions of less than 5 mK per day. However, because the 
underground laboratory is only at a depth of 5 m, there 
is still a peak to peak temperature variation of about 1 
degree per year, accounting for the change of seasons. 
Changes in the atmospheric pressure also affect the di- 
mensions of the ring laser structure by changing the com- 
pression of the Zerodur block and cannot be neglected. 
Hence G is kept in a pressure stabilized enclosure. A feed- 
back system based on the determination of the current 
value of the optical frequency of the lasing mode of one 
sense of propagation allows for active control of the pres- 
sure inside the steel vessel such that an overall geometric 
scale factor stability of better than 10"^*^ is routinely ob- 
tained. At the same time the design of the instrument is 
made as symmetric as possible. So changes in area and 
perimeter are compensated with a corresponding change 
in wavelength as long as no shear forces are present and 
the longitudinal mode index stays the same. 

A typical eight day long measurement sequence of rota- 
tion rate data from the G ring laser is shown in Fig. [8j In 
order to demonstrate the obtained sensor sensitivity we 
have subtracted the mean Earth rotation rate from the 
gyroscope data. The y-axis gives the measured variation 
of the rate of rotation, while the x-axis shows the time 
expressed in the form of the modified Julian date. Each 
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FIG. 8: Approximately eight days of raw G data taken with 
30 minutes of integration time. One can clearly see the con- 
tributions from diurnal polar motion, solid Earth tides and 
local tilt. 



data point was taken by integrating over 30 minutes of 
measurement data. There are several distinct signal con- 
tributions in the data, which come from known geophys- 
ical effects. The most prominent signal is caused by di- 
urnal polar motion (25j . The polar motion data is super- 
imposed by a tilt signal caused by the semi-diurnal and 
diurnal tides of the solid Earth, distorting the otherwise 
sinusoidal diurnal frequencies slightly. At the Geodetic 
Observatory in Wettzell the tilt effects of the solid Earth 
tides can be as large as 40 nrad in amplitude. In Fig. [8] 
the diurnal signal is dominated by the polar motion |36| 7 
Less evident in Fig.|8]are the effects from local tilt, which 
contains periodic signals of tidal origin as well as non- 
periodic signals. The latter are non-periodic and usually 
change slowly over the run of several days. High res- 
olution tiltmeters inside the pressure stabilizing vessel 
of the G ring laser keep track of these local effects and 
the data is corrected for gravitational attraction (atmo- 
sphere, sun and moon)l [25j. Large non-periodic local 
tilts occur most prominently after abundant rainfall, in- 
dicating hydrological interactions with the rock and soil 
beneath the ring laser monument. Fig. [9] shows the east 
component of three tiltmeters installed i) on a gravime- 
ter pillar at the surface, ii) in 6 to depth, and iii) in 30 
TO depth. 

While the tilmeter in 30m depths clearly shows the 
periodic signal of the solid earth tides, the tilt record 
of the instruments near to the surface is dominated 
by large non-periodic signals hydrological, thermoelastic 
and barometric origin. Several investigations have shown 
that the site and the installation depths of tiltmeters has 
a major impact on environmental noise mainly coming 
from hydrology [?T], [35], [33], [32] has shown that even 
in 100 TO depths effects caused by hydrological changes 
are detectable, but strongly reduced in comparison to a 
50 TO deep installation. First investigations related to to- 
pographic and temperature induced effects were carried 
out by [32] and [30]- Detailed investigations using the 
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FIG. 9: Measurement of local tilts as a function of depth in 
the Earth. 
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finite-element method have shown that these effects can 
amount to more than 10 nrad (|37j. |38j). while the dis- 
tance between the source and the location of observation 
can be several hundred meters. Additionally, recent work 
using the G ring laser data reveals that effects caused by 
wind friction at the Earth surface yields to high frequency 
rotations of large amplitudes. 

The large seasonal temperature effect on the G ring 
laser as well as the substantial local tilt signals and the 
rather high ambient noise level of our near soil surface 
structures give reasonable hope of much better perfor- 
mances of a ring laser installation in a deep underground 
laboratory such as the Gran Sasso Laboratories. 

For the detection of fundamental physics signals one 
has to remove all known perturbation signals of the Earth 
from the ring laser time-series. Furthermore we have ap- 
plied 2 hours of averaging of the data in order to reduce 
the effect from short period perturbations. Fig. [T0| shows 
an example. In Fig.[TT]we show the current sensitivity ex- 
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FIG. 10: The rotation rate of the Earth measured with the 
G ring laser as a function of time. Averaging over 2 hours was 
applied to a corrected dataset, where all known geophysical 
signals have been removed. 

pressed in term of Allan deviation of the G, the expected 
sensitivity of each ring laser at Gran Sasso Laboratories 
and the relevant geophysical signal. 

In order to reduce the local orientation uncertainties, 
which remain after local tilts measured with the high 



FIG. 11: Resolution and stability of G, compared with Earth 
signals 




FIG. 12: Six rings arranged on the faces of a Cube, using the 
GEOSENSOR design, which has been successfully used so far 
for middle size rings, as our prototype G-Pisa 



resolution tiltmeters have been removed, averaging as in- 
dicated above was applied to a series of 30 days of data 
collection, including the period shown in Fig. [Sj It can 
be expected that a similar data set from the Gran Sasso 
laboratory would become substantially smoother, since 
most of the perturbations, caused by ambient atmosphere 
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- topsoil interaction still contained in the data of Fig 
would no longer exist in the deep underground facil 
ity. Changing hydrologic conditions presumably causing 
small local rotation and temperature variations, atmo- 
spheric pressure and wind loading are among the sources 
for the systematic signatures in the residual data. 



V. CONFIGURATION OF A TRI-AXIAL 
DETECTOR 

From now on, we will restrict our analysis to 24 m 
perimeter rings, arranged in two configurations that are 
of some experimental interest, i.e. 6 ring lasers rigidly 
mounted on the faces of a cube, as shown in Fig. [12] 
and 3 ring lasers oriented along the edges of an octahe- 
dron, see Fig. [l3] The cubic configuration requires 24 
mirrors forming 6 independent rings and the extension 
of the GEOSENSOR design is straightforward (see sub- 
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FIG. 13: Three rings are formed using 6 mirrors located on 
the vertices of an octahedron 



section VC ); while the octahedral configuration require 
6 mirrors only to form 3 orthogonal rings. By itself the 
configuration which uses a cube is redundant, each ring 
has a parallel companion, which can be used for the study 
of systematics. For the octahedron configuration the im- 
plementation of the GEOSENSOR design needs further 
development. Redundancy can be easily obtained con- 
structing a second octahedron with planes parallel to the 
other one. The two structures should be built very close 
to each other, in order to keep as much as possible the 
whole apparatus compact; in this way 6 rings are avail- 
able, analogously to the cube configuration, see Fig. [14] 
This configuration has the advantage that there are con- 




FIG. 14: Six rings, two by two parallel, with mirrors on the 
vertices of two octahedron, constructed very close one to the 
other in order to reduce the dimension of the apparatus. 

straints in the relative angle between rings, since each 
mirror is in common between two rings, and three linear 
Fabry-Perot cavities are available using the three diago- 
nals of the rings. Those linear cavities have the capability 
of monitoring the relative angles between different rings, 
and as well the length of each diagonal. 



A. Ring-laser sensitivity 

The rotation sensitivity ctq, for noise fluctuations 
which are dominated by laser shot noise over an inte- 
gration time T, reads 



hf 



cP 



AAQ V WT 



(28) 



where Q is the quality factor of the optical cavity, / = 
c/A is the laser frequency, h is the Plank constant and 
W is the power of the laser [27] . The limiting sensitivity 
can be conveniently calculated scaling the parameters of 
the WettzeU "G" ring laser 




In order to obtain in few weeks a 10% accuracy level 
in the measurement of the relativistic effective rotation 
rates, we must achieve the sensitivity goal of do = 
7 X 10"^'* rad/s (or equivalently a rotation noise level 
20 prad/ sec/ Hz^/^ at a frequency of 1 day~^). From Eq. 
p9| we have that a system of 6 rings with P = 24 m, 
g = 3 X 1012 and W = 200 nW can fulfill this require- 
ment. 



B. Expected performances of not optimally 
oriented rings 

We assume that the ring lasers are identical in the 
sense described in Sect. IIII El and that the dihedral an- 
gles arccos(MQ • up) are measured better than one part 
in 10^° in order to estimate fl independently from the 
reference frame. Note that only the stability of dihedral 
angles can be monitored by means of the Earth signal it- 
self only for short times (few days), while their measure- 
ments and controls must be performed independently in 
the laboratory. For instance, assuming that the scale fac- 
tors are controlled to the 10~^'^ accuracy, the responses 
of two parallel rings is statistically different from noise 
when their parallelism is modified. 

From an experimental point of view, to arrange in the 
Cartesian planes several rings and keep the configuration 
stable over the integration time T ~ 1 day is a demanding 
task. However, we can relax such a demanding require- 
ment by means of data analysis procedures that account 
for slightly non-orthogonal dihedral angles. 

For instance, we can use the measured dihedral angles 
to estimate directly Q. In fact, we can substitute the 



quadratic combination of ring laser responses in Eq. ( 17 ) 
with the equivalent bilinear combination 



M 

E 



(30) 
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where Qap are the elements of the M x M matrix Q = 
N{NN'^)-'^N'^. The statistics of Q is no longer non- 
central x^; however, we can easily compute (see App. |b] 
for details) its mean 

< Q >= \n\^ + Maf, (31) 

and variance 

a'g = 24 J2 + J2 Ql^ulul (32) 

In the limit of high SNR, fluctuations of Q tend to be 
Gaussian distributed, and so we recover the results in Eq. 
( [22| ) for the overall sensitivity of the system. If we start 
with dihedral angle close to 7r/2 (say 1 part in 10^), then 
sensitivity loss is very small since it is of the same order. 




FIG. 15: Drawing of G-Pisa, based on the GEOSENSOR 
design 



C. Guidelines of the Experimental Apparatus 

The best performing ring, so far, is G which is a four 
mirrors ring. This is one of the reasons why the present 
scheme uses a square ring geometry. In principle a trian- 
gular ring, with 3 mirrors could be preferable since the 
three mirrors are always inside a plane, and the losses will 
be minimized as well, reducing the number of mirrors. 
It could be advantageous in principle, but a triangular 
ring is less sensitive. For instance, let us compare the 
performance of two rings inscribed in a circle of radius 
r; for a regular polygon with different number of sides 
N, the area is ^ = ^ sin(^) and the perimeter is 
P = 27V r sin(^); it is straight forward to demonstrate 
that the triangular ring has 0.7 times the signal than the 
square one, which is equivalent to say that the triangular 
ring needs 2 times more time to reach the same level of 
accuracy as the square one. 

The ring-laser response is proportional to the Scale Fac- 
tor " S" . For a perfect square ring this proportionality 
factor is equivalent to N the number of wavelength inside 
the ring: when the length of the ring changes, because of 
a change in the temperature, the laser changes its wave- 
length in order to keep N constant. This is true as long 
as the perimeter change is below a wavelength, 632 nm 
in our case, and in this conditions the gain factor of the 
instrument guarantees a very high accuracy of the mea- 
surement. For example: if the laboratory has ST — 1° 
degree temperature excursion, the ring perimeter is 36 
m, in order to guarantee the operation of the ring-laser 
with a fixed number of wavelengths N, it is necessary 
to realize the whole apparatus using materials with tem- 
perature expansion coefficient of the order of 10^® K^^. 
This is the concept used for G in Wettzell: a structure 
realised with material as Zerodur, with a design which 
can be defined monolithic, i.e. relative motions of the 
mirrors are not allowed. G has a very high stability, but 
is rather expensive, and not very fiexiblc with regards to 
changing the mirrors and align the laser cavity. More- 
over the extension of this design to a large array of rings 



seems rather difficult. Later on, a more flexible and less 
expensive design has been realized, called GEOSENSOR, 
which so far has been employed especially for smaller size 
rings. This design allows a very good relative alignment 
of the mirrors, it is relatively easy to change mirrors and 
tools to move each mirrors along different degrees of free- 
dom have been implemented. So far this kind of instru- 
ments have been done in steel. Fig. [15] shows a drawing 
of G-Pisa, our prototype. The optical cavity vacuum 
chamber has a stainless steel modular structure: 4 tow- 
ers, located at the corners of the square and containing 
the mirrors holders inside, are connected by pipes, in or- 
der to form a ring vacuum chamber with a total volume 
of about 5 • 10"'^ m^. The mirrors are rigidly fixed to the 
tower. The cavity alignment can be adjusted by moving 
the towers with respect to the slab through a lever system 
that allows 2 degrees of freedom of movements. No win- 
dow delimits the active region and the vacuum chamber 
is entirely filled with a mixture of He and a 50% isotopic 
mixture of ^°Ne and ^^Ne. The total pressure of the gas 
mixture is set to 560 Pa with a partial pressure of Neon of 
20 Pa. The active region is a plasma produced in a cap- 
illary pyrex tube inserted at the middle of one of the ring 
sides by a radio frequency capacitively coupled discharge. 
In a non monolithic device, temperature changes could 
interrupt the continuous operation, and the perimeter is 
actively controlled by acting on the mirrors and using as 
reference a stabilized laser; very highly stabilised lasers 
are commercially available, for instance wavelength sta- 
bilization at the level of 2.5 x 10^^ using iodine line can 
be obtained. G-Pisa is kept in continuous laser opera- 
tion trough a perimeter stabilization servo system which 
acts along the diagonal direction, for two opposite placed 
mirrors, through piezoelectric actuators [52] . 
The GEOSENSOR design has other advantages as well: 
the mirrors are under vacuum and are not affected by the 
outside pressure changes, they can be very easily aligned 
and the cost is pretty much reduced compared with the 
monolithic design. The experience of G-Pisa has shown 
so far that it can work with different orientations. In fact 
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FIG. 16: Mace, for a rectangular ring, with sides 6 x 6.6 m, in 
function of a misalignment of one of the four mirrors with re- 
spect to the ideal position, the perimeter control acts on four 
mirrors, maximum thermal excursion of 1 degree and the sup- 
port of the GEOSENSOR has thermal expansion coefficient 
of 7 X 10"'' K'^ (granite) 



G-Pisa has worked both horizontally and vertically ori- 
ented. It is in steel, inside the thermally stabilized room 
in the central area of Virgo, in order to improve thermal 
stability, it has been mounted on top of a granite table 
(thermal expansion coefficient about 5 x 10"^ m/mK). 
To find the guidelines of the mechanical project, we have 
used a simple program which consists in considering the 
ring as four points (the light spot on the mirrors) which 
can be moved from the ideal position, both inside the 
plane or outside the plane. The model takes into account 
thermal expansion and the perimeter is kept constant by 
acting diagonally on pairs of mirrors; the use of 2 mirrors 
or 4 mirrors for the feedback correction have been inves- 
tigated; the thermal excursion is considered of 1° degree. 
The scale factor S in presence of misalignments is com- 
pared with 5*0 ( scale factor at the optimal configuration) ; 
this comparison is expressed as M^cc = ^"sa^ i which gives 
the accuracy limit induced by misalignments. The re- 
quired level of accuracy of 1 part of 10^° is Mace = 10~^°. 
Fig. [16] shows Mace for a rectangular ring, with sides 6 
m and 6.6 to, in function of a misalignment of one of the 
four mirrors. 

Fig. [l6] clearly shows that the Gain Factor changes a lot 
with small change of mirrors positions. The situation 
strongly improves by considering a perfect square ring. 
In fact, for a closed figure with a fixed number of sides, 
the area over perimeter ratio has a maximum when the 
polygon is regular one, as for example a 'perfect' square 
ring. Fig.[l7|and[l8|show Mace with 100 ^ra construction 
precision and two possible choices of the thermal expan- 
sion coefficient. For instance, let us assume that each 
mirror position is in the ideal position within a quan- 
tity i5 which depends on the precision of the construc- 
tion. Fig. 19 shows Mace when 3 out of the 4 mirrors are 



FIG. 17: Mace in function of a misalignment of one of the 
four mirrors with respect to the ideal position, the perimeter 
control acts on four (thick line) or two (dashed line) mirrors, 
maximum thermal excursion of 1 degree and the material has 
7 X 10"'^ m/mK (granite) 

Thermal expansion 10 times less than granite 
1 degree themperature change 
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FIG. 18: Same as Fig. |17| but with expansion coefficient 10 
times lower 



positioned with an error, 10000 points have been evalu- 
ated pseudo-randomly distributed between ±50/xto along 
each coordinate. In summary, if the thermal excursion is 
1 K, the position of the mirrors is an ideal square within 
±50/im, the support has a thermal expansion coefficient 
below 7 X lO^^iiT"^, Mace remains in the range necessary 
for the needed accuracy using four or two mirrors active 
control of the perimeter. 

Misalignments which bring the light spots outside the 
plane of the ring do not have appreciable effect on the 
gain factor, but they change the orientation of area vec- 
tor Mq,; in this case the effect for Mace depends on the 
relative angle between the ring and the Earth rotational 
axis. Figs. 20 and 21 shows how the accuracy changes 



for two different ring orientations: parallel to the axis of 
the Earth and at 45° degrees respectively. The first is 
almost insensitive, while the other is sensitive to nano- 
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FIG. 19: Histogram of Mace when the position of 3 mirrors 
are within ±50 /im close to the ideal position, 10000 points 
have been evaluated by randomly extracting the position error 
(±50/im). The thermal expansion coefficient is 7 x 10~^A"~^, 
thermal excursion 1 K, top histogram shows the case with 4 
mirrors control of the perimeter, bottom curve with 2 mirrors 
control. 



Area vector 45 degrees with respect to the Earth axis 
granite thermal expansion 

Relative Accuracy Error 




FIG. 21: Accuracy change in percentage for vertical misalign- 
ments and area vector close to 45 degrees with respect to the 
Earth rotational axis. The ring geometry is not perfect in the 
plane, there is a misalignment of lOO/im, a maximum temper- 
ature change of 1 K, and the thermal expansion coefficient is 
7 X 10"'^m/m/_K' 



Area vector parallel to the Earth axis 
granite thermal expansion 

Relative Accuracy Error 
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FIG. 20: Accuracy change in percentage for vertical misalign- 
ments and area vector close to the parallel alignment to the 
Earth rotational axis. The ring geometry is not perfect in the 
plane, there is a misalignment of lOOfim, a maximum temper- 
ature change of 1 K, and the thermal expansion coefficient is 
7 X W'^m/m/K 



Relative accuracy of the Earth angular rotation 
in function of the inclination angle 

Relative Accuracy Error 




FIG. 22: Relative limit of the accuracy in the measurement 
of the Earth angular rotation induced by a nm change in 
the position of one of the mirrors with respect to its original 
position, in function of angle with the Earth rotation axis. 
The area vector of the ring lays in the meridian plane. The 
accuracy loss is zero when the Earth axis and the area vector 
are parallel, and is very high in the orthogonal alignment 



metric misalignments. 

Fig. [22] shows Mace for a nm vertical misalignment of 
one of the rings in function of the angle with respect to 
the Earth rotational axis. 

In summary: The Gain Factor of each ring can be kept 
constant at the level of 1 part in 10^°, if the positions of 
the mirrors are constructed and kept within +/ — 50 /im 
error close to the ideal square ring; the relative position 
between mirrors can be rigidly constrained with gran- 
ite, super-invar or similar low thermal expansion coefh- 
cient spacers, it is preferable to use all the four mirrors 
for the perimeter active stabilization, but two mirrors 



control could be acceptable as well if the structure has 
thermal coefficient better than granite. It is necessary 
to constantly monitor the relative angle between rings 
with nrad precision (only the relative alignment mat- 
ters). This can be accomplished looking at the modal 
structure of the FP cavities formed along the diagonals. 
Moreover the orientation of each ring with respect the 
Earth rotation axis should be such to avoid alignment 
too sensitive to the relative angle (relative angle with the 
Earth rotation axis below 60°). Using the Earth angular 
velocity rotation, which is perfectly stable for few days, 
the whole apparatus can be calibrated at the beginning; 
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the relative angle, or the area of each ring could be not 
perfectly planar or exactly 90", but it is important to 
monitor the geometry of the structure during the whole 
measurement time (years). The mirror holders plays an 
important role, it can be advantageous to build them in 
Zerodur or similar material, in order to avoid displace- 
ments out of the plane. The mirrors holder should be 
designed in order to provide the tools to align the cav- 
ities; in principle each mirrors should have 5 degrees of 
freedom: three translations and two tilts, the rotation 
around the axis orthogonal to the mirror itself does not 
play a role; but since the mirrors are spherical only three 
motions are fundamental: we may have one translation 
along the diagonal and two mirrors tilts or three transla- 
tions. 

Let us consider now an octahedral geometry, containing 
the three rings. 

Fig.[22]shows that the relative angle between the differ- 
ent rings must be monitored at the level of nrad, and that 
it should be avoided to put one of the rings with an angle 
larger than 60" with respect the Earth rotation axis. We 
have done the exercise to fit the octahedron, with rings of 
24 m perimeter, inside the node B of LNGS, considering 
that this node is 8 m tall, and imposing the constrains 
discussed in Fig. [22j The exercise is done with the oc- 
tahedron since it needs more space. Considering that 
the latitude of LNGS is 42° 27" N two configurations are 
given: to have the octahedron straight up (8.4 m tall) or 
laying on one side, one ring is respectively horizontally 
or vertically oriented and the other two symmetrically 
positioned with respect to the meridian plane. Let us 
consider the maximum size 9 m: 8.48 m, the diagonal 
of the octahedron, plus 0.6 m necessary to hold mirrors 
and optics in general necessary for the read out. This 
octahedron can be contained inside each of the big halls 
of LNGS, in both orientation, but inside node B, which 
is the most isolated room of LNGS, see Figs. [23] and [25] 
the only possibility is the shorter configuration, with the 
longer side parallel to the floor. 

So, the octahedron, with rings of 24 m perimeter, can 
be contained inside node B, where the ceiling is 8 m tall, 
while for node C the structure should be scaled, probably 
no more than 20 m perimeter can be contained inside 
node C, since the ceiling there is 6 m tall. Fig. [24] shows 
the octahedron inside node B. 




FIG. 23: Plan of LNGS laboratory close to node B 




VI. DIAGNOSTICS OF DIHEDRAL ANGLES 
AND SCALE FACTORS 

To reduce the influence of systematics in long-term 
measurements, the control of the geometrical stability of 
ring laser system is of paramount importance. In partic- 
ular, it is crucial to monitor the deviations from planarity 
of each ring laser and their mutual orientations. 

A square ring consists of four spherical mirrors with the 
same curvature radius R, placed at the corners. Square 
geometry guarantees that opposite mirrors are parallel so 



FIG. 24; the octahedron inside node B 



that they form two extra linear Fabry-Perot cavities (see 
Fig. ~ 



26). As a consequence, each square ring is made 



of three optical resonators: the ring itself and two linear 
ones oriented along the diagonals. These latter can be 
used to monitor the geometrical stability of the whole 
ring system. Deviations from a square geometry result 
in tilting and/or displacements of the diagonal vectors, 
which in turn change the cavity eigenmodes. 

A linear symmetric FP cavity with spherical mirrors in 
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|C£m|"' / |Coo|^ 


l = Q 


1= 1 


1 = 2 


m = 


1 


5.16 X lo'e^ 


1.33 X lo'-'e;^ 


m — 1 


5.16 X lo'e^ 


2.66 X lo'^'e^ey 




m = 2 


1.33 X lO'^'e* 







TABLE I: Power coupled to the first cavity higher modes 
(£ + m = 0,1,2) as a fraction of tlie external laser power for 
Qc ~ Qx = Qy and X = Y = 0. The value are obtained for 
a ratio between the cavity length and the mirror radius of 
curvature of 1.5 



FIG. 25: The ring laser system inside node B of LNGS, side 
view showing that passage between the two entrances 



Zm = i|c2 = {L being the square ring arm) and 

centers on the z-axis, supports the Gaussian modes 



1 



^/2a; ' 
{zY^' \ Wc {z) 



H„ 



Wc (z) 



xexp 



2;2 _|_ y2 / 2^ 

-ik — ikz + i (£ + m + 1) arctan — 

2qc(z) \ b 



(33) 



where qdz) ^ z ^ ib = - z and b = 

y/d {2R - d) are the complex curvature of the Gaussian 
beam and the confocal parameter, respectively; here the 
curvature radius Rc {z) and the spot-size Wc {z) read 



Rc {z) = 

Wc (z) = 



d^ - 2dR - 4z2 
Iz 

A 4z2 + 2dR - 
TT 2^d{2R-d) 



The eigenmodes Ei^m {x, y, z) form a complete set which 
can be used for representing a generic field confined be- 
tween the two generally misaligned mirrors of the cavity 

E{x,y,z) ^'^Cl^rnEl^rn{x,y,z) , 



e,m 



where 



dx J dyEin (x, y) Ei^^a [x, y, zm) 



and Ein {x, y) is the beam illuminating the input mirror 
Ml . If we suppose the mirror tilted by Qx and Qy and 
displaced by X and Y with respect to cavity axis z, we 
have 



— ik 



E,r,{x,y)o^e [^'"^(^m) ^""i'^') 



r+e. 



As an example, the relative intensities \Ctm\ /I Cool 
for the first modes t + m — 0, 1, 2 and X = Y = Q are 
reported in Tab. |l) 

It is clear that the cavity axes misalignment can be 
detected by looking at the intensity pattern of the beam 



transmitted though the output mirror M2. A modal de- 
composition of such a pattern gives a suitable set of coef- 
ficients \Cim\ which can be used for estimating the posi- 
tion and angular misalignment of the cavity with respect 
to the reference beam E^ {x,y). 

Supposing that at the begining it — 0) the cavity ex- 
ternal laser is perfectly aligned to a symmetric cavity 
(if the two mirror show equal transmittivity then the 
cavity transmission is 1) so that all the incoming power 
Pin is coupled to the TEMOO mode. The measurement 
procedure we have devised is a tunable laser, showing a 
linewidth narrower than the cavity linewidth, tuned over 
a cavity FSR in a time interval At so that each mode is 
spanned in a time t = where is the cavity finesse. 
The number of photons in the (.m mode are given by (we 
are now assuming a rectangular line shape instead of a 
Lorentzian profile) 



nir, 



hv 



This number of photons must be higher than the noise 
equivalent number of photons hitting the detector in 
the same time interval. The noise equivalent power in 
W/ \fTTz^ is connected to the equivalent number of pho- 
tons by 

NEP^fB 
nNEP = — ; T 



hv 



where t] and B are the quantum efficiency and the detec- 
tion bandwidth, respectively. 

To overcome the photon noise, we have to satisfy the 
inequality 



riNEP njsiEP |C( 



00 1 



, < 1 

In particular, looking at Cqi coefficient we have 

1 nepVb 

"™ ^ 5.16 x 10762"^ ' 
further, by assuming rtoo — j^''' we obtain 

NEPs/B 



5.16 x 10^927; 
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For typical silicon detectors NEP - lO'^'^W/v^, B - 
l{fHz, 7] - 0.9, so that 

P,n > 2.2 X IQ-^^Q-^ . 

For a tilt sensitivity of 8 ^ 10^^ the power required at 
the input is 

P,„ > 220 mW . 



10"" rad/sec] here we have used GM = 3.986004418 x 
10^4 m3/s^ R = 6.378137 x 10^ m and Ts = 
86164.0989 s. Assuming one year of data taking with 
the same ring laser parameters used for the simulations 
in Sect. HID we have that the standard deviation of ft'g 

and Q'j. is an^g — 0.03, and therefore upper limits of 
some interest can be put on ai and 7 at the Gran Sasso 
colatitude ~ tt/A. 




FIG. 26: In a square ring configuration passive Fabry-Perot 
cavities are formed along the square diagonals (dashed line in 
the sketch). In the case of an octahedron each of these three 
passive cavities is shared by two rings. 



VII. MORE ABOUT THE IMPORTANCE OF 
THIS MEASUREMENT 

A. Post Newtonian Parameters 

The proposed experimental apparatus is well suited for 
performing optical test of metric theory of gravitation. 
We start from the statement that the vector fl' should be 
entirely contained in the meridian plane if the preferred 
frames effect, determined by W (see Eq. (A9)), can be 
neglected. Indeed the currently available best estimates 
P5] suggest that this effect is about 2 orders of magnitude 
smaller than the geodetic and Lense-Thirring contribu- 
tions. As a consequence, we expect that the measured 
components of f2' outside the meridian plane should be 
compatible with noise. In this case, our results could 
be used to obtain new constraints, independent from the 
available ones, on the preferred frames parameters. In 
addition, we can write the PPN parameters ai and 7 as 
a function of the Ur and ue components of CI' 

ai = (-4fi^ CSC d~8^^ sec 9 

j-i = (n'gcsce-n[.sece/2j ~2 , (34) 

where fl'^g = fl'^g/w and w is the very precisely mea- 
sured constant w = 2ttGM/{c'^RTs) ~ 5.0747798 x 



B. Interdisciplinary: Geodesy and Geophysics 

Earth rotation rate and the orientation of the rota- 
tional axis of the Earth in space are the linking quantities 
between the terrestrial (ITRF) and the celestial (ICRF) 
reference frames. Currently a set of quasars, forming an 
external set of markers, provide the only way of deter- 
mining the rotational velocity and the variations of the 
orientation of the rotational axis of the Earth with suf- 
ficient accuracy. As already mentioned, 10 fis for the 
measurement of length of day (LOD) and 0.1 mas for 
the pole position are routinely achieved by a network of 
VLBI radio telescopes as one of the services (lERS) of 
the International Association of Geodesy (lAG). The op- 
eration of such a network requires expensive equipment 
and a lot of maintenance effort. Huge amounts of data 
are recorded in each measurement session, which require 
physical transport over large distances for the correlation 
in the analysis centers. Data latency and the fact that 
there is no continuous measurement coverage are suggest- 
ing the investigation of alternative methods for the pre- 
cise estimation of Earth rotation. Furthermore it is desir- 
able to develop an independent measurement technique, 
in order to identify intra-technique biases if they exist. 
Ring lasers are possible candidates for such an alterna- 
tive measurement technique. They measure the earth 
rotation locally and within much shorter time intervals. 
Such gyros are widely used in aircraft navigation and can 
measure rotations absolute, i.e. independent of an exter- 
nal reference frame. Therefore also local contributions to 
earth rotation are contained in the measurements. The 
effects of earth tides, strain, crust deformation, seismic 
events, polar motion are contained in the ring laser mea- 
surements due to their contribution to earth rotation or 
due to variations in the orientation of the respective ring 
laser. However, the demands on such instruments are ex- 
tremely high and cannot be met by existing commercial 
devices. They can be summarized as: 

• sensitivity to rotation 0.01 prad/s at about 1 hour 
of integration 

• sensor stability of 1 part in 10^'^ over several month 
to years (Chandler Wobble) 

• resolution in sensor orientation w 1 nrad. This 
corresponds to polar motion of around 1cm at the 
pole. 
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This means that a reasonable improvement in sensor sen- 
sitivity and stabihty is stiU required in order to make 
ring lasers viable tools to be applied to space geodesy. 
The design of the G ring laser is one way of approaching 
these demands and it not too far away from reaching this 
goal [45j . Operating several such ring laser gyroscopes in 
geophysical independent regions simultaneously offers a 
unique possibility to distinguish global from local (nion- 
umentation related) signal contributions through their 
independent data streams. 



Discussions and Conclusions 

The feasibility of the experiment for the measurement 
of relativistic effective rotation rates appears to rest only 
on a tri-axial dynamical sensor of local rotation of enough 
sensitivity. Despite the fact that large ring lasers as G are 
very stable platforms and with the provision of tight feed- 
back systems to stabilize the scale factor (cold cavity, as 
well as the active cavity), currently ring laser gyroscopes 
are not able to determine the DC part of the Earth rota- 
tion rate with a sensitivity compatible with the require- 
ments for detection of the Lense-Thirring effect. While 
the contribution of the varying Earth rotation itself pre- 
sumably can be removed with sufficient accuracy from 
the C04 series of VLBI measurements, there remains the 
problem of determining the actual null-shift offsets from 
the laser functions in the ring laser gyroscope. Since 
the gravito-magnetic effect is small and constant, a good 
discrimination against laser biases, such as for example 
'Fresnel drag' inside the laser cavity must be achieved. 
Therefore it will be advantageous to add one or several 
ring laser cavities in addition to the triad structure for 
sufficient redundancy. We also intend to operate at least 
the G ring laser structure in parallel to the here proposed 
structure in order to discriminate local perturbation sig- 
nals from regional and global ones. A second large ring 
laser located at the Cashmere facility in Christchurch, 
New Zealand, will be used in the data analysis process, 
provided it can be run with sufficient resolution and sta- 
bility. 



Appendix A: Ring Laser Measurements in the 
Laboratory Frame 

In this Appendix we evaluate the response to the grav- 
itational field of a ring laser in an Earth bound labora- 
tory and, to know the space-time metric in the laboratory 
frame we shall use the construction of the "proper refer- 
ence frame" as described in Ref . |71 |3S] ■ 

As we discussed in Section III, a ring laser converts a 
time difference into a frequency difference (see e.g. Eq. 
([S])). It is possible to show that (see e.g. [55]) in a sta- 
tionary metric in the form|59| 17^1/ = gf^^^x^) an observer 
at rest at = Xg measures the proper-time difference 
(5t = T+ — r_ between the right handed beam propaga- 



tion time (t_|-) and the left handed one (t_): 

St = -2./5oo(4) f —ds' = -2./goo(4) f H ■ ds, 
* Js 900 * Js 

(Al) 

where S is the spatial trajectory of the beams, whose 
tangent vector is ds, and we set Hi 



goi 



In order to evaluate the proper-time difference (All, 



we need to know the space-time metric in our laboratory, 
that is to say the gravitational field nearby the world- 
line of the observer which performs measurements with 
the ring laser. To this end, we consider an observer in 
arbitrary motion in a given background space-time, and 
write the corresponding local metric in a neighborhood 
of its world-line (see e.g. [7]) 



f(0)(0) 

5(o)W 



= l + 2A-x + 0{x'^), 
- r!(,)(,.)a:W+0(a;2), 
= ?7W(j) + 0{x^). 



(A2) 
(A3) 
(A4) 



It is worth pointing out that the Eqs. (A2)-(A4) hold 



only near the world-line of the observer, where quadratic 
displacements terms are negligible. Here we suppose that 
the observer carries an orthonormal tetrad (parentheses 



whose four-vector e 



(0) 



co- 



refer to tetrad indices) 6(0), 
incides with his four- velocity U, while the four- vectors 
e(i) define the basis of the spatial vectors in the tan- 
gent space along its world-line. By construction we have 



where 



is the Minkowski ten- 



e(Q)e(/3) = 

sor. The metric components (|A2p-(|A4p are expressed 



in coordinates that are associated to the given tetrad, 
namely the space coordinates cc''^ and the observer's 
proper time x^°\ In the above equations, A is the spa- 
tial projection of the observer's four-acceleration, while 
the tensor r2(i)(fc) is related to the parallel transport of 
the basis four-vectors along the observer's world-line: 



"(a) = -t^(/3)^ (a)- -^^^ parucuiar, n 
the tetrad would be Fe rmi- Walke r transported. Let us 
remark that the metric ( A2 )-( [A4| ) is Minkowskian along 
the observer's world-line (?*•' =0); it is everywhere flat 
iff >t = 0, i.e. the observer is in geodesic motion and the 
tetrad is non rotating (i.e. it does not rotate with respect 
to an inertial-guidance gyroscope). In the latter case, the 
first corrections to the flat space-time metric are 0{x^) 

In order to explicitly write the local metric, which 
through its gravito-magnetic (goi) and gravito-electric 
(goo) components enables us to evaluate the proper-time 
difference ( Al ), we must choose a suitable tetrad by tak- 



In particular, if were zero. 



ing into account the motion of the Earth-bound labora- 
tory in the background space-time metric. To this end, 
we consider the following PPN background metric which 
describes the gravitational field of the rotating Earth (see 
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e.g. 



(1 - 2U{R))dT^ - (1 + 2-fU{R)) S^jdX'dX^ + 



(1 + 7 + q;i/4) 

i?3 



(JeAi?),-aiC/(i?)W, 



dX'dT, 
(A5) 



where —U{R) is the Newtonian potential, is the an- 
gular momentum of the Earth, Wi is the velocity of the 
reference frame in which the Earth is at rest with re- 
spect to mean rest-frame of the Universe; 7 and ai are 
post-Newtonian parameters that measure, respectively, 
the effect of spatial curvature and the effect of preferred 
frames. The background metric (A5) is referred to an 



Earth Fixed Inertial (ECI) frame, where Cartesian geo- 
centric coordinates are used, such that R is the position 
vector and R = ^/T^ = ^/ X"^ + + Z'^ . Then, 
we choose a laboratory tetrad which is related to the 
background coordinate basis of ( A5 ) by a pure Lorentz 



boost, together with a re-normalization of the basis vec- 
tors: in other words the local laboratory axes have the 
same orientations as those in the background ECI frame, 
and they could be physically realized by three orthonor- 
mal telescopes, always pointing toward the same distant 
stars. 

In this case, one can show that the gravito-magnetic 
contribution in the local metric reads [71 HHl HH ISO] 
^{i)(k)x'^''^ — — (rJ' A cc)(.j, where the total relativistic 
contribution 17' is the sum of four terms, with the di- 
mensions of angular rotation rates 



defined by 
J7t = 



-(l + 7)VC/(i?) AV, 
l + 7 + ai/4 /Je _ 3J^R 
2 \R^ i?5 

ai^VU{R) A W, 
1 , dV 



R 



(A6) 

(A7) 
, (A8) 

(A9) 
(AlO) 



The vector fJ' represents the precession rate that an 
inertial-guidance gyroscope, co-moving with the labora- 
tory, would have with respect to the ideal laboratory spa- 
tial axes (see e.g. [71 135]) which are always oriented as 
those of the ECI frame; if the spin vector of the gyroscope 
is S, its precession is hence defined by 



dS 

~dt 



n' AS 



(All) 



Differently speaking, we may say that the local spatial 
basis vectors are not Fermi- Walker transported along the 
laboratory world-line. In particular the total precession 
rate is made of four contributions: i) the geodetic or de 
Sitter precession Q,q is due to the motion of the labora- 
tory in the curved space-time around the Earth; ii) the 



Lense-Thirring precession fls is due to the angular mo- 
mentum of the Earth; iii) Jln/ is due to the preferred 
frames effect; and iv) the Thomas precession fix is re- 
lated to the angular defect due to the Lorentz boost. 

It is worth noticing that for a laboratory bounded to 
the Earth 



A^'^-VUiRl 



(A12) 



and the acceleration A. can not be eliminated. However, 
for a geodetic motion (e.g. a free fall satellite) ^ = 
and by substituting in Eqs. (A7| and (AlO), we obtain 



and 



Og = - 



1 



+ 7 VU{R)AV, 



1 



Aav = o. 



(A13) 



(A14) 



Eq. (A13) gives the geodetic precession for a free fall 



gyroscope, while Thomas precession (A14| is zero. 
All terms in (|A7|)-( AlO ) must be evaluated along the 



laboratory world- line (hence, they are constant in the lo- 
cal frame) , whose position and velocity in the background 
frame are R and V, respectively. However, if we consider 
an actual laboratory fixed on the Earth surface, the spa- 
tial axes of the corresponding tetrad rotate with respect 
to the coordinate basis of the metric ( A5 1 , and we must 



take into account in the gravito-magnetic term ( A3 1 the 



contribution of the additional rotation vector fl^ , which 
corresponds to the Earth rotation rate, as measured in 
the local frame [60]. 

As a consequence, it is possible to show that, up to lin- 
ear displacements from the world-line, the relevant local 
gravito-magnetic potential turns out to be 



9(o){^) = (^^ Aa;)(^) 



(A15) 



where fl = — — Jl', while the gravito-electric ff(o)(o) 
one remains the same. 

Now, we are able to evaluate the proper-time difference 



St 



'2ygao{xl) f H -ds 



(A16) 



Without loss of generality, we suppose that the observer 
is at rest i n th e origin of the coordinates, so that, ac- 
cording to (A2|, 500(2^0) = 1- As a consequence, we have 



6t 



{n A x) 
s (1 + 2^- a;) 



• ds. 



(A17) 



Now, on taking into account the expression of acceler- 
ation of the laboratory frame (A12) and evaluating the 



magnitude of the various terms, the leading contribution 
to (A17l con be written, applying Stokes theorem 



(5r = 



-2 [ [V A (17 A x)] ■ dA, (A18) 

J A 
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where A = Aun is the area enclosed by the beams and 
oriented according to its normal vector m„. On evaluat- 
ing the curl, taking into account that Jl is constant, we 
eventually obtain 



St = -4 n-dA = -An ■ A. 



(A19) 



On substituting fl — —fl^ — fl' in (A19), we see that 



the proper-time delay can be written in the form 
St = 4^35 • A + 4rj' • A, 



(A20) 



where 4170 • A is the purely kinematic Sagnac term, due 
to the rotation of the Earth, while Afl' ■ A is the gravi- 



tational correction due to the contributions (A7)-(A10) 



According to Section III A from Eq. (A19I, it is then 



possible to write the ring laser equation in the form 



Sf 



AA 
XP 



n. 



(A21) 



To further clarify Eqs. (A7)-(A10) it is useful to use an 



orthonormal spherical basis Ur,u^, in the ECI frame, 
such that the d — tt/2 plane coincides with the equa- 
torial plane. As a consequence, the position vector of 
the laboratory with respect to the center of the Earth 
is i? = Ruj. and the kinematic constraint V = A R 
holds, i.e. V = flmRamOui^. 



Thus, the components of Jl' in physical units read 

fir - 



c/R 

1 + 7 + q;i/4 G 



2 

ai GM 

1 
2^? 



c'^RS 

Ur A W, 



[Je -3(Je ■Ur)Ur], 



^1%R^ sin^ iJfie, 



(A22) 
(A23) 
(A24) 
(A25) 



r 



Moreover, we assume the general relativistic values of the 
PPN parameters, 7 = 1, ai = 0, and use for the Newto- 
nian potential of the Earth its monopole approximation. 



i.e. U{R) = GM/R. Thus, the components ( [A22| )-( [A23| ) 

read 



GM 

-2—— sin -fffl^Ui}, 
c^R 

G 



Je - 3 (J© • Ur) Ur 



0, 



2c2 



^%R^ sin^i9fJe> 



(A26) 

(A27) 
(A28) 
(A29) 



and, to leading order, the total rotation rate which 
enters the Eq. (|A20|) is 



O = Of, 



GM 

2—— sin i?ri®Me 



G 



[J® ~ 3 (J© • Ur) Ur] 

(A30) 



If we denote by a the angle between the radial direction Ur and the normal vector it„ , on setting m„ = cos aUr 
sinaug in (A20), and using (A3G|, we may express the proper-time delay in the form 



St = 



4A 



Gla 



a) — 2^^r — il© sin ^? sin a + " ^ fj© (2 cos cos a + sin sin a) 
c^R c'R-^ 



17© cos ( 

where we have written Jm = ImClm, in term of the /©, the moment of inertia of the Earth. 



(A31) 



r 



Appendix B: Probability Distribution of Quadratic if a; is a multivariate Gaussian random vector with mean 
Forms s and covariance matrix S, the mean and the variance 



The statistics of quadratic forms of Gaussian random 
vectors x are well known in the literature. In particular. 
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of a quadratic form Q — x^Qx are given by 

<g> = <x'^Qx>=Tr{Q'E) + s^Qs 
<tI = < {x^Qxf > - < Q >2 



(Bl) 



where Q is a square symmetric matrix, ^ and Tr are the 
transpose and trace operators, respectively. The statis- 
tics of Q in general is not known, unless QS is an idem- 
potent matrix |5Ij. In the case were x represents the 
response of ring lasers in a regular polyhedral configura- 
tion Q = I, with no common noise source and the same 
sensitivity S = cr^/, where I is the identity matrix, the 
above formulas greatly simplifies 



< g > 

2 



Ma^ + E 
2Ma^ + AEa^ 



(B2) 
(B3) 



where E — s*s — ||s|p is the signal energy. In this case 
also the statistics of Q readily follows. In fact, starting 
from its definition we have 

Pq{Q) = j P{x)5{Q^x'^x)dx (B4) 



where P{x) = exp[(a; - s)^ {x - s) / {2cj'^)]/ {2TTa'^)^/'^ 
is the Gaussian probability density of one sample of the 
random vector x. We can use the integral representation 
of the Dirac's 5-function 



r+oo 

5{Q-x'x) = / e- 



(B5) 



and write 



Pq{Q) 



+00 



dwe*"'^ 



1 



exp 



(27ra2)M/2 
1 



iuj x^x — -T—^ix — s)^ Ix — s) 



Ax 



By re-arranging the exponent, the last integral can be re- 
cast as a M-dimensional Gaussian integral and calculated 
explicitly 



(27ra2)W2 



exp 



1 

2(j2 



T 
X X 



T 
tS X 



E 



da; 



(1 - 2ia;r72)M/2 



(B6) 



where in the last expression one recognizes the moment 
generating functions of a non-central distributions 
with AI degrees of freedom and non-centrality param- 
eter E. The probability density function of Q can be 
found using the tables of Fourier Transform pairs 



written using their Fourier transforms, 



+00 



duj e" 



= -exp[-{Q + E)/{2a^)] 



exp[iu}cr^E/{l — 2iuja^) 
(1 - 2iwa2)W2 

M-2 

Q 



where Ik{x) are the modified Bessel functions of order k. 



Appendix C: Probability Distribution of Projectors 

The norm of complementary projection operators P 
and Q acting on Gaussian random vectors x are de- 
scribed by remarkably simple statistics. In fact, starting 
from the definition of iJp = ||Pa;||2 a.nd Eq = WQxW^ we 
have that the joint probability density P{Ep, Eq) reads 

P{Ep,Eq)^ J P{x)5{Ep-x^Px)S{EQ~x^Qx)dx 

(CI) 

where P{x) is the probability density of one sample of 
the random vector x. The two Dirac (5-functions can be 



P{Ep,Eq) ^ J p(a;)eK^+"-E«-^''("-P+'''3)^] dudv dx 

(C2) 

where the integrals in du and dv are performed along 
the imaginary axis (i.e. u = iuji and v = 1002 are purely 
imaginary complex numbers). Now assume the noise is 
Gaussian distributed, uncorrelated between different de- 
lM/2-ii\/ QE/a'^)),ectoTS and with identical variance cr^ in every detector, 
namely 



Pix) = 



1 



(2^2)^ 



exp 



1 

2^ 



{x - sf{x ~ s) (C3) 



where s = (Jl • i6i, . . . r2 • u^^^) is the rotation signal in 
vectorial form. Then, 



P{Ep,Eq) = [-) ' j Gicp[-a{x^sf{x^s)\ X 

X exp [— a;"'"(uP -f wQ)a;] da; • 

X e"-^^e"-^« dudu, (C4) 

where a = 1/2(T^. Writing a; as s + e and switching the 
integration variable to e yields 



22 



j exp [-'nFial + uP + vQ)n-2ri^[uP + vQ)s\ deexp [-s^ {uP + vQ)s\ q^Ep^vEq 



The integration in An can be done by noting that it is a 
standard Af -dimensional Gaussian integral with the lin- 
ear term, and in general, for any M x M symmetric ma- 
trix A and M-vector b, 



j exp(- 



A n + n) dn 



M/2 



ydet(A) 



exp 



In our case. 



A al + uP + vQ 
b = 2{uP + vQ)s. 



b'^A-^ 



(C6) 



(C7) 



Now we exploit the properties of P and Q. Using their 
complementarity, we can write 



A = {a + u)P +{a + v)Q 



(C8) 



and from the fact that they are orthogonal and idempo- 
tent we also have 



A-' 

hence 



-s^Ps 



a + V J 



(C9) 



(CIO) 



Furthermore, as P and Q are projection matrices, their 
eigenvalues are {0, 1} with multiplicities respectively 
{M - 2, 2} for P and {2, M - 2} for Q. Then, writing 
A in diagonal form is trivial and leads to 

dct{A) = {a + uf{a + v)^'-^, (Cll) 

determinants being independent from the basis. By us- 
ing 



CIO and Cll in C6 one can see that the Gaussian 



integral splits into the product of factors involving either 
u or V. By further substituting in |C6[ the remaining in- 
tegrals separate and the probability density remarkably 
factorizes as 



(C12) 



with 

PiEp) = 
P{Eq) = 



1 



+ 20-2^ 

1 

1 + 2a'^v 



exp 



-So U 



1 - 

— -1 



2a\ 



du 



exp 



2a^v 



dv 
(C13) 



and Sb 



Ps. 



s^Qs. The transformed functions 



is the moment generating functions of two non-central 
distributions, with 2 and M—2 degrees of freedom respec- 
tively, and whose non-centrality parameters are P s 
and s^Q s respectively. Thus, 



P{Ep) 
P{Eq) 



2 exp 



exp 



Ep 



2^2 

Eq + 



2^2 



Eq 

Sr, 



(C14) 



where J„(a;) is the modified Bessel function of the first 
kind. Some interesting conclusions can be drawn about 
the virtual channels Ep and Eg, which make them in- 
teresting for the identification of meridian plane and the 
estimate of fl. 



1. Ep is distributed as a non-central with 2 degrees 
of freedom and non-centrality parameter equal to 
s^P s, i.e. the magnitude of the signal projection 



in the P subspace. 



2. Eq is distributed as a non-central with M — 2 
degrees of freedom and non-centrality parameter 
equal to s^Q s, i.e. the magnitude of the signal 
projection in the Q subspace. 



3. Ep and Eq are statistically independent processes. 



4. In the limit of high SNR, Ep and Eq are Gaus- 
sian distributed with means < Ep >— s^P s, 
< Eq >= s^Q s and variances a'^^ = Aaf^s'^P s, 
a% — {M — 2)(TqS^Q s, respectively. 
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